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ABSTRACT 

Observed Hot Jupiters exhibit a wide range of physical properties. For a given 
mass, many planets have inflated radii, while others are surprisingly compact 
and may harbor large central cores. Motivated by the observational sample, 
this paper considers possible effects from collisions of smaller rocky planets with 
gas giant planets. In this scenario, the Jovian planets migrate first and enter 
into ~ 4 day orbits, whereas rocky planets (mass nip = 0.1 — 20M^) migrate 
later and then encounter the gaseous giants. Previous work indicates that the 
collision rates are high for such systems. This paper calculates the trajectories of 
incoming rocky planets as they orbit within the gaseous planets and are subjected 
to gravitational, frictional, and tidal forces. These collisions always increase the 
metallicity of the Jovian planets. If the incoming rocky bodies survive tidal 
destruction and reach the central regions, they provide a means of producing large 
planetary cores. Both the added metallicity and larger cores act to decrease the 
radii of the gas giants at fixed mass. The energy released during these collisions 
provides the Jovian planet with an additional heat source; here we determine the 
radial layers where kinetic energy of the colliding body is dissipated, including 
the energy remaining upon impact with the existing core. This process could 
have long-term effects if the colliding body deposits significant energy deep in 
the interior, in regions of high opacity. Both Hot Jupiters and newly formed 
gas giants have inflated radii, large enough to allow incoming rocky planets to 
survive tidal disruption, enhance the central core mass, and deposit significant 
energy (in contrast, denser giant planets with the mass and radius of Jupiter are 
expected to tidally destroy incoming rocky bodies). 



Subject headings: planets and satellites: composition — planets and satellites: 
dynamical evolution and stability — planets and satellites: formation 
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Introduction 



Since the first discovery of an exoplanet orbiting a main- sequence star ( iMayor fc Queloz 



19951 ). tlie field of extrasolar planets lias undergone rapid progress. More than 750 planets 
have been discovered beyond our Solar Sy stem, with over 1000 more possible objects awaiting 
further measurements for confirmation (IBorucki et al.ll201ll ). Although the observational 
sample is far from complete, it is clear that the Galaxy hosts a diverse population of planets, 
with a range of masses spanning over three orders of magnitude. More than 150 of the ~ 750 
confirmed planets are gas giants with semi-major axes less than 0.1 AU. These planets are 
thought to reach such extre me environments after forming further out in their circumstellar 



disks and migrating inward ( jPapaloizou fc Tergueml 120061). Up to ^ 5% of all solar systern s 



are believed to harbor "Hot Jupiters" (IBramich &: Hornd 120061 : Ivan Saders fc Gaudil 1201 ll ) 
so that a billion such planets could reside in the Galaxy. 



The observed Hot Jupit ers have a wide range of radii for a given mass ( iBaraffe et al. 



2010l : iBurrows &: Ortorul2009l ). Many planets are larger than predicted, but some are surpris- 
ingly small. The reasons for these anomalous radii remain uncertain, and the construction of 
a definitive m ass-radius theory for Hot Jupiters has become a central problem in exoplanet 



research (e.g.. iBodenheimer et al.ll2003uLaughlin et al.ll201ll ). In genera l terms, the inflated 



planets can be unde rstood if they contain additional energy s ources jBodenheimer et aL 
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Gu et al 



2010 



20041). e.g., heat gener ated by Ohmic dissipation (IBatygin fc Stevenson 



Perna et al.boid iBatvgin et al.ll201l[) , conversion of wind k inetic energy into heat ([Rauscher &: Menoul 
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i 



Guillot fc Showman 2002), and/or tidal dissipation ( Arras &: SocratesI 2010). In 



con- 



trast, giant planets with deflated radii raust h ave high metallicities ( iMiller &: Fortneyll201l[ ) 
and possibly large cores ^Burrows et al. 2007). The problem of these anomalous radii thus 
splits into two sub-problems: [a] the source(s) of additional heating, and [b] the ways in 
which planets can accumulate metallicity and large cores. Identifying possible long-term 
heat sources and examining the ways in which planets can increase their core masses will 
help explain the observed mass-radius relation and increase our understanding of giant planet 
formation, structure, and evolution. 

The goal of this paper is relatively modest: We consider the effects of collisions of 
smaller rocky bodies with gaseous giant planets, and explore how the collisions affect metal- 
licity, core size, and perhaps the energy budget. During the early epochs of solar system 
formation, a large number of small terrestrial bodies are present; previous work has shown 
that they can rnigrate inward and become captured by a Hot Jupiter with high probability 
( iKetchum et al.ll2011bl ). After capture, a rocky planet orbiting within the interior of a Hot 
Jupiter loses energy due to dissipative forces and eventually either colhdes with the core or 
becomes tidally disrupted. These impacts provide a natural mechanism for increasing the 
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planetary metallicity, and could help explain some of the inferred heavy cores and metal-rich 
compositions of the Hot Jupiters. In addition, a substantial amount of energy is released 
during these collisions. This energy could inflate a giant planet over a longer period of time, 
provided that the energy is deposited deep within the planet interior where the opacity is 
high. On the other hand, if the kinetic energy from the collision is dissipated in the outer 
planetary layers, it would provide only a short-term source of energy. These impacts thus 
have the potential to help explain, at least in part, the observed radius anomalies. 

Given the strong observational evidence for planets with large cores, the process of 
core enhancement is of particular interest. A specific observational example is provided 



by the planet HP 149026 
0.73 ± 0.03i ? 7 JSato et al. 



calculations (IFortney et al. 



wit h a total mass Mp = 0.36 ± 0.03Mj and a radius Rp 



2006 



2005 ) . These observati onal constraints, along with theoretical 
Ikoma et al.ll2006[ ). indicate that this planet contains ~ 80 



of heavy elements. A significant fraction of these heavy elements could reside in the core. 
Although current observations do not fully constrain the core mass, it is expected that HD 
149026b and similar planets have larger cores than predicted by theories of planet formation. 
According to the core accretion model (jPoUack et al.l 119961 ). giant planet cores are formed 
from rocky planetesimals in the circumstellar disk, and begin to accrete gas after reaching a 
critical mass around lOM^. Sufficiently large planets experience runaway gas accretion and 
clear gaps in their disks, so that the accretion of additional rocky material is suppressed, 
making it difficult to build up core masses as large as 80M®. The planet HD 149026b, and 
others, motivates the possible role of collisions in shaping giant planet cores. A gaseous 
planet with an initially small core mass could accrete one or more rocky planets early in its 
lifetime and thereby accumulate a heavier core. Collisions with superearths {nip ^ lOM®) 
are probably the best way to explain the most massive cores, since the accumulation of a 
large core from low-mass planets would require a high encounter rate. In addition, small 
planets are unlikely to survive tidal forces over the entire trajectory and reach the core. 
However, the accretion of smaller planets could affect the structure of the giant planets by 
increasing the metallicity of the envelope. 

In this paper we model the coUisional dynamics between Jovian planets and smaller 
rocky p lanets. This problem has beeri explored previously in the context of Jupiter and 
Saturn ( lAnic et al.l 120071 : iLi et al.ll2010[ ) using smooth particle hydrodynamical (SPH) sim- 
ulations. These authors have investigated both the effects on the overall structure of the 
giant planet and survival of the colliding bodies. The conditions required for rocky planets 
to survive are particularly important. The previous results indicate that low-mass plan- 
ets {nip ^ IM®) are often destroyed prior to reaching the core, whereas larger planets 



[nip ~ 



IOM0) usually survive the entire trajectory. This paper extends this previous work, 
but adopts a complementary approach. SPH (or equivalent) methods provide a detailed de- 
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scription of the collisions, but are computationally expensive and allow relatively few cases 
to be considered. Here, instead, we treat the collision as a one-body problem, where the in- 
coming rocky body orbits within the giant planet and is subject to gravitational, frictional, 
and tidal forces. This approach allows for a much wider survey of parameter space, which, 
in spite of the simplifying approximations, increases our physical understanding of these 
collisions. More specifically, this paper investigates the accretion of rocky bodies spanning 
two orders of magnitude in mass (mp = 0.1 — 20M^), within gas giant planets of varying 
properties. Here we focus on Jovian planets with mass Mp = IMj a nd radii in the rang e 
Rp = l — 2Rj, where the latter are appropriate for young Hot Jupiters (IBurrows et al.lll997l ). 

This paper is organized as follows: In Section [2l we construct polytropic models for 
gas giant planets both with and without central cores. Using this model in Section [31 we 
calculate the trajectories of incoming rocky bodies within the gaseous giant planets. These 
results determine whether the rocky body is captured or passes through the giant planet (as 
a function of initial speed and impact parameter). These trajectories also specify the radial 
locations where energy is dissipated and the kinetic energy remaining upon impact with the 
existing core. In Section HJ we calculate the tidal forces exerted on an incoming rocky planet 
by its host Jovian planet; such tidal effects are a function of radial distance and determine 
whether or not the rocky planet can survive. We conclude, in Section [5l with a summary of 
our results and a discussion of their implications. 



2. Planetary Structure Models 

Before considering planetary collisions, we develop a structure model for the Jovian 
(target) planet. Here we assume that the Jovian planet consists of a rocky core of constant 
density surrounded by a gaseous envelope. The planet is assumed to be in hydrostatic 
equilibrium, and obeys the equations of hydrostatic equilibrium and mass conservation, 

dP GM{r)p{r) ^ dM . ^ ^ 

— = ^ ^ and — — = 47rr>(r) . (1) 

dr dr 

After specifying the equation of state P{p), these equations can be combined into a single 
differential equation, which can be solved for the density profile. The conditions within 
Jovian planets require complicated equations of state, which can vary from planet to planet. 
Since the goal of this paper is to develop a relatively simple model to survey the parameter 
space for collisions, we adopt a polytropic equation of state of the form 

P = Kp^ where 1 = 1 + -, (2) 

n 
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and where n is the polytropic index. Here we consider a range of polytropic indices n = 1 — 2, 
but detailed models sugge st Jovian plane ts are characterized by indices near the lower end 



[n ~ 



1) of this range (e.g., lBurrowsll201ll ). As shown below, the mean density of the planet 



is one of the most important physical variables in the problem; for polytropic models with 
fixed mass and radius, the mean density does not differ greatly from that of more detailed 
models. Further, the system to system variations in the structure of observed exoplanets is 
larger than the differences resulting from this polytropic approximation. 



Following the formulation of IChandrasekhan (119391 ). we define dimensionless quantities 
and /(^) through the relations 



^ = - and p = Pc/"(0 > where 

JrC 



KT 



(3) 



With these definitions, the resulting differential equation for /(^) becomes the Lane-Emden 
equation 

1 ^ / _^f\ 

(4) 



1 d 



with the boundary conditions 



dC 



0, 



(5) 



where C,c denotes the radius of the planetary core. Planets without central cores correspond 
to the limit ^ 0. The first zero of the function /(^), denoted here at ^o, marks the outer 
boundary (the surface) of the planet. 



Although polyt ropic stellar structure models have been used extensively (j Chandr asekhar 



19391 : |Phillipslll994l ). this application differs from most previous cases in two respects: First, 
we include the possibility of a constant density core, so that the (dimensionless) core radius 
C,c provides an additional parameter. Second, the expected values of the polytropic index are 
closer to n = 1 for giant planets, instead of the more familiar value n = 3/2 for degenerate 
stars and convective stars, or n = 3 for massive stars dominated by radiation pressure 



flPhillipsI 119941 ). 



In terms of the dimensionless variables, the mass conservation equation is given by 

dp 



2 rn 



(6) 



which is related to the physical enclosed mass profile M(r) via 



M(r)=4rfpcM0- 



(7) 
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The total dimensionless mass /^o of the planet is thus 

/io= r^nOd^, (8) 

^0 

and the physical mass and radius are given by 

Mp = AnR^p^fio and Rp = R^q. (9) 

Although the planetary structure is hydrostatic and does not require a time variable, 
orbital motion is time dependent. Here we define a dimensionless time variable 

T = - where to = (47rGpJ~^/^ . (10) 

^0 

The time scale to is the free-fall time appropriate for the density at the planet center, and 
provides a benchmark dynamical time scale for the orbit calculations. As an order of mag- 
nitude estimate, the central densities of Jovian planets pc ~ 10 g cm~^, which corresponds 
to a time scale to ~ 345 sec ~ 6 min. 

Numerical integration of equation (jll) yields a dimensionless density profile, which can 
be converted to physical units by specifying the central density. The surface of the planet 
corresponds to the radius where the density and pressure vanish. Solutions for /(^) 
are shown in Figure [1] for a polytropic index n = 1 and dimensionless core radii = — 1 . 
Variation of the parameters and n lead to a range of density profiles, and upon specification 
of the physical planet mass and radius, different central densities and pressures, core masses, 
and core radii. This model thus allows us to explore a variety of planet structures. For this 
paper, we vary the core radius over a wide range, with a focus on = — 1, and vary the 
polytropic index over the range n = 1 — 2. For planets with Jovian masses and radii, we 
obtain central densities pc = 2 — 15 g cm~'^, leadi ng to central pressur es of order 100 Mbar, 



in agreement with more complicated models (e.g., iMilitzer et al.ll2008l ). For a dimensionless 
core radius = 1, the core mass (in physical units) thus has values mc ~ 20 — 30 (for 
varying n). 

In general, equation (|4]) must be numerically integrated. For the particular case where 
the polytropic index n = 1, however, the Lane-Emden equation has an analytic solution 
of the form 

/«) = ^S|1 + B2|i. (11) 

For a coreless polytrope, the first term must vanish to avoid divergence at the origin [A = 0), 
and the solution becomes 

/«) = . (12) 
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0.2 0.4 0.6 0.8 1 

Fig. 1. — Density profiles for n = 1 polytropic models with a collection of core masses. Each 
curve corresponds to a different core radius from = (bottom) to = 1 (top) . For the 
largest value shown, = 1, the core makes up ~ 6.6% of the mass (e.g., the core would have 
mass rric ~ 21Me for a giant planet mass Mp — IMj). 
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Note that the first zero occurs at ,^0 = tt. For finite core radii {C,c 7^ 0), the general solution 
takes the form 

/(O = (^ccos^c - sin^c) ^^Y^ + (^csin^c + cos^J (13) 

Note that for this general case, the first zero ^0 niust be determined by solving a transcen- 
dental equation. However, analytic expressions for the mass profiles /i(^) can be found by 
substituting these solutions for /(^) into equation ([6]). For the coreless polytrope, we find 

/t^(0 = cose + sine (14) 
For the general case, = C,^/3 for ^ < ^c, and 

= l^c' + i^c - cos (e - ^c) + (1 + U) Sin (e - ^c) (15) 

for C, > Cc- Although the coreless polytrope represents a relatively good model for plan- 
ets with small cores like Jupiter, it is less accurate for the exoplanets with massive cores. 
Nonetheless, it provides a reasonable approximation and the analytic forms for the density 
and mass profiles are convenient for later calculations. 

For a coreless polytrope of fixed index n, the model provides a two parameter family 
of solutions, where the parameters are the central density pc and the coefficient K of the 
equation of state (note that K depends on the physics of the planetary interior and can take a 
range of values). For any planet with observed mass Mp and radius Rp, we can (in principle) 
adjust pc and K to provide the correct properties. This procedure results in a prediction 
for the central density pc and the central pressure Pc = Kp^. As a consistency check, we 
fit the model to a planet with the mass and radius of Jupiter, calculate the central pressure 
and density, and then compare with current estimates. For example, a coreless model with 
n = 1 produces a central density pc ~ 4 g cm~^ and a corresponding central pressure 
pressure Pc ~ 40 Mbar. The actual central properties of Jupiter are somewhat uncertain 
because the interior cannot be observed directly and the conditions are difficult to recreate 
experimentally. However, the central pressure of Jupite r is thought to lie in the range Pc 



= 10 - 100 Mbar, with estimates around Pc ~ 70 Mbar (jSaumon fc Guillotll2004j ). in rough 
agreement with our predicted model values. More detailed knowledge of these conditions 
requires a reliable equation of state for hydrogen at high pressures, especially in the transition 
to metallic hydrogen. This poses significant experimental and theoretical cha.llenge s, and 



many uncertainties remain regarding the interiors of gas giants (iHubbard et al.ll2002l ) 
The central pressure of the planet must obey the relation 



p. = :f^ = KpI. (16) 
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The first expression follows from the condition of hydrostatic equilibrium (and dimensional 
analysis) whereas the second follows from the polytropic equation of state. The dimensionless 
parameter T is expected to be of order unity and can be written in the form 

T= ^ (17) 

For an n = 1 polytrope with no core, = = /iq, so that T = vr^/S. As the core size 
increases, the factor J-" decreases, so that J-" ^ 0.242 for a core with C,c = 2. For larger values 
of polytropic index n, the factor J-" is larger for planets with no core, but has nearly the same 
values for planets with large central cores. 



3. Orbital Trajectories 

In this section we solve the equations of motion for a rocky planet as it falls into a 
gaseous giant planet. In this setting, we assume that the Hot Jupiter has completed its 
inward migration and entered into a stable orbit close to its star. Rocky planets migrate 



in later and have a high probability of coll iding with the Hot Ju 



orbital eccentricity is not overly damped (iKetchum et al.l l2011bl ). At the start of these 



3iter, provided that their 



integrations, the rocky planet already lies well within the Hill sphere of the Jovian planet 
so that the gravitational field of the star can be neglected. The system thus consists of a 
Jovian planet with a mass Mp = IMj and a smaller rocky planet with mass rrip <^ Mp. To 
first approximation, the system can be treated as a one-body problem, with the origin of 
the coordinate system placed at the center of the Jovian planet. As the rocky planet orbits 
within the gaseous envelope of the Jovian planet, it experiences the usual gravitational force 
and a frictional (ram pressure) force, so that the equation of motion is given by 

GM{r) ^ p{r)v^j: 



. , , p 
r 1 



where Sp ~ nrp^ is the effective cross sectional area of the rocky planet (for determining 
ram pressure forces). Using definitions from equation (|3]), we obtain 



R ^ 



-f - aPv'v , (19) 



where v is now the dimensionless speed (in units of R/tQ), and where time derivatives are 
performed with respect to the dimensionless time r (from equation [10]). We have also 
defined a dimensionless friction coefficient 

iL (Ml) (1} 



"-T^^l^ll^l- (20) 
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Note that the coefficient a depends on the masses and radii of both planets, as well as 
the polytropic index and the core radius. Here we consider rocky planets with masses 
rrip = 0.1 — 20M^. For a typical mean density (p) = 5.5 g cm~^ for the rocky planet, and 
for Mp = Mj and Rp = Rj, the corresponding values of the coefficient a lie in the range 
a ~ 0.5 — 5. Note that Jovian planets with inflated radii have lower values of a and provide 
less friction. 

The relevant parameter space thus consists of five variables: the impact parameter xq 
of the collision, initial velocity Vq, friction coefficient a, the polytropic index n, and the 
dimensionless core radius ^c- Once these parameters are specified, the equation of motion 
f|T9|) can be numerically integrated. 



3.1. Rocky Planets Orbiting within Jovian Planets 

To begin, we consider the simplest case of a coreless Jovian planet (^c = 0) and adopt 
a polytropic index n = 1. By fixing the polytropic index and core radius, the number of 
free parameters is reduced to three. The starting conditions are defined as follows: The 
calculations are started with the incoming rocky planets located at the outer edge of the 
Jovian planet, traveling in the —y direction. In dimensionless coordinates, the rocky planet 
has initial position {xo,yo), where Xq is identified as the impact parameter, so that the 
starting coordinates obey the constraint 

4 + yl = ^0 • (21) 

Here we explore the full range of impact parameters < xq < ^o- For the sake of definiteness, 
we consider incoming rocky planets with masses rrip = 1 and 10 M^. If the Jovian planet has 
the mass and radius of Jupiter, these values corres pond to dimensionless friction coefficients 



in the range a = 1 - 3. Dynamical simulations (iKetchum et al.ll2011bl ) indicate that the 
expected impact velocities lie in the range f o = 40 — 150 km s~^. This study considers a 
wider range of impact velocities where f o = 1 — 10^ km s~^. 

Examples of typical trajectories are shown in Figure [2] for a dimensionless impact pa- 
rameter Xq = 1.5, an initial velocity vq = 40 km s~^, and rocky planet masses rup = 1 and 
10 (friction parameters a = 1 — 3). As the rocky planets spiral inward within the Jo- 
vian planets, the orbits of larger rocky planets remain relatively circular, while the orbits 
of smaller rocky planets become highly elliptical. In general, the bodies continuously spiral 
inward and approach the origin asymptotically. 



There are two possible outcomes of these orbits: Either the rocky planet is captured 
and spirals inward, or it passes through the gaseous planet and retains enough kinetic energy 



.4L . . . J 

-4 -2 2 4 
X 



-0.2 0.2 

X 




Fig. 2. — Orbits within coreless planets with polytropic index n = 1. The top panels show 
results for incoming rocky planets with masses rUp = IMq; and bottom panels show results 
for rrip = lOM®. The panels on the left show a full planetary view, whereas the panels on 
the right show the orbits on a smaller spatial scale near the center. Both collisions use an 
impact parameter xq = 1.5 and an impact velocity vq = 40 km s^^. 
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Fig. 3. — Penetration velocity as a function of impact parameter. The velocity fpen is the 
impact velocity (in km s~^) needed for the rocky planet to pass through the gaseous giant 
planet, emerge from the far side, and remain uncaptured. The impact parameter is given 
as a fraction of the planetary radius. Here the giant planet has mass IMj, radius IRj, and 
n = 3/2. The solid curve shows the numerical results; the dashed curve shows the analytic 
approximation derived in Appendix |Al The approximation works well over the full range of 
impact parameters, with the largest discrepancy for head-on collisions {xq = 0). 
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to escape. For an incoming rocky planet with given impact parameter xq, we can determine 
the initial velocity necessary for the rocky planet to avoid capture. More specifically, the 
critical condition for penetration requires that the rocky planet initially has enough kinetic 
energy so that it passes through the gaseous planet and leaves its surface traveling outward 
at the escape speed. This penetration speed is denoted here as fpen and is a function of xq- 
For a given impact speed, incoming rocky planets will be captured for all impact parameters 
Xq < x*, where the maximum impact parameter defines the effective capture cross section 
a{vo) = Tcxl. 

Here we calculate the required penetration speeds numerically using the above formula- 
tion. We also derive an analytic estimate of the penetration speed Vpen in Appendix |X1 Figure 
in] shows the numerically determined penetration speeds along with the analytic approxima- 
tion. The best agreement between the numerical and analytic results occurs for large impact 
parameters because the assumption of constant density (see Appendix [A]) is more accurate 
far from the center of the planet; even for small impact parameters, however, the results 
agree within a factor of two. These calculations show that enormous starting velocities are 
required for penetration when the impact parameter is near zero, i.e., fpen ~ 10^ km s~^. 
For somewhat larger impact parameters Xq ~ ^o/2, penetration still requires fpen ~ 1000 
km s~^. As a result, the cross section for capture is large. For the expected impact speeds 
vq = 40 — 150 km s~^, incoming rocky bodies will penetrate the Jovian planet only for impact 
parameters comparable to the radius of the planet, roughly xq > 0.8.^o = x*- For high-speed 
rocky planets with vq > 100 km s~^, the effective target area (capture cross section) of the 
Jovian planet is thus reduced by ~ 30%. 

We now consider the general case where rocky planets collide with Jovian planets con- 
taining central cores. For the sake of definiteness, we begin by choosing a standard (dimen- 
sionless) core radius = 1-0. Depending on the polytropic index, this core radius spans 
~ 20 — 30% of the giant planet radius, and encloses a mass of 20 — 30M^. Here we use three 
polytropic indices {n = 1, 3/2, and 2), four terrestrial planet masses {nip = 0.1, 1, 10, and 
20 M®), and the full range of impact parameters xq = — ^c, where the initial velocities 
span the range Vq = 20 - 150 km s~^. The resulting trajectories, including profiles of energy 
dissipation and the velocities at the core, are nearly independent of the polytropic index. 
Given that similar results were obtained for all three choices {n = 1, 3/2, and 2), only the 
results for the n = 1 polytrope are shown here. Figure H] shows the trajectories for rocky 
planets with initial speeds of 80 km s~^ and with varying impact parameter. 

We note that the presence of a core has little effect on whether or not an incoming rocky 
planet will penetrate the Jovian planet and escape capture. As a result, the results depicted 
in Figure |3] remain valid — only incoming planets with large impact parameters {xq > O.Si^o) 
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Fig. 4. — Trajectories for different impact parameters. Each curve shows the orbit of an 
incoming rocky body for a given impact parameter. The outer dashed circle marks the 
location of the outer surface of the Jovian planet; the inner dashed circle marks the location 
of its core, which has dimensionless radius = 1. The impact speeds, defined at the surface 
of the Jovian planet, are taken to be 80 km s^^. 
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or incredibly large initial speeds {vq > 1000 km s ^) are able to escape. 

3.2. Energy Dissipation 

During each incoming trajectory, energy is dissipated through the action of frictional 
forces. Figures E] and E] show the total energy (kinetic and potential) of the rocky planet as a 
function of radial distance, measured from the center of the giant planet. These figures show 
results for target giant planets with the mass and radius of Jupiter, and for incoming rocky 
planets with mp = IM^ (top panels) and nip = IOMq (bottom panels). Figure E] shows 
energy profiles for varying initial impact speeds, whereas Figure E] shows energy profiles for 
varying impact parameters. In all cases shown, the majority of the energy is dissipated in 
the outer layers {r/Rp > 0.5) of the planet. 

As shown in Figure El for sufficiently large impact parameters the energy fraction E/Eq 
takes on multiple values for a given radius; on the other hand, the radius is a single-valued 
function of energy. This complication arises for orbits that reach an inner turning point 
and reverse course before finally spiraling inward (e.g., see the orbit with the largest impact 
parameter in Figure H]). For rocky planets with larger masses and/or giant (target) planets 
with lower density, the orbits display additional turning points and the energy dissipation 
curves show even more structure (these cases are not shown in the figures). 

When the rocky planet strikes the core, its motion is arrested and any remaining ki- 
netic energy is converted into heat. The energy remaining upon impact with the core thus 
determines the extent of possible core deformation as well as the amount of energy available 
for long-term heating. This fraction of the initial energy remaining depends on the mass of 
the rocky planet and the structure of the Jovian planet. Rocky planets with larger masses 
retain a larger fraction of their initial energy (as well as a larger absolute value of energy). 
For target giant planets with the properties of Jupiter (see Figures [5] and E]) , 1 planets 
lose most of their initial energy, whereas 10 planets retain up to 30%. The structure 
of the Jovian planet, in particular its mean density, plays the most important role. When 
the radius of the giant planet is doubled (which may be appropriate for young planets), 
the 10 planets can retain the majority of their initial energy (not shown). An analytic 
desription of the energy dissipation process is given in Appendix [Bl 

The kinetic energy of the rocky planet at the moment of impact with the core is shown 
in Figure [7] as a function of both the initial velocity and impact parameter. The impact 
velocities vary over the range 10 — 75 km s~^, leading to impact energies in the range 
Iq38 _ ]^g42 gj^gg^ where these values are strongly dependent upon the rocky planet mass. 
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Fig. 5. — Fraction of the total energy of the rocky planet as a function of radial distance 
for varying initial velocities. Masses for the rocky planet are rup — 1 (top panel) and 
rup — lOM® (bottom panel). The Jovian planet has polytropic index n — 1, mass Mp — IMj, 
and radius Rp — IRj. The dashed hues indicate the core radius. 
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Fig. 6. — Fraction of the total energy of the rocky planet as a function of radial distance 
for varying impact parameters. Masses for the rocky planet are nip — 1 M® (top panel) 
and nip — lOM® (bottom panel). The Jovian planet has polytropic index n — 1, mass 
Mp — IMj, and radius Rp — IRj. The dashed lines indicate the core radius. 
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Planets less massive than rrip = IM® show little spread in their final energy (velocity) for 
the entire set of impact parameters and initial velocities. This result can be understood in 
terms of the nature of the frictional force {Ff ~ pv"^) and the size of the friction coefficient 
a. Planets with the highest velocities experience the highest levels of damping, and their 
speeds are quickly diminished. Furthermore, the impact parameter determines the average 
density encountered by the orbiting planets. Although the average density is lower for planets 
incident at large impact parameters, these planets must travel farther before being accreted 
onto the core, and are thus subjected to damping forces acting over larger distances. These 
effects are the most pronounced for low-mass planets, since they experience the greatest 
frictional forces relative to gravitational forces. Planets with masses larger than rrip ~ IM^ 
are subjected to relatively less friction and are more sensitive to variations in the initial 
velocity and impact parameter. Planets incident at small impact parameters (roughly, less 
than the core radius) collide directly with the core and travel shorter distances. These planets 
thus experience less damping and retain a significant amount of energy upon impact with the 
core. Incoming planets with sufficiently large impact parameters do not collide directly with 
the core, but rather gradually spiral inward and often enter into highly elliptical orbits. The 
orbital path length thus increases substantially for impact parameters larger than a critical 
value, and the amount of energy deposited at the core decreases accordingly. For the regime 
of parameter space considered here, this critical impact parameter is roughly xq* ~ OARp. 
Finally, we note that the curves in Figure [7] do not extend beyond xq ^ 0.9Rp] for collisions 
with larger impact parameter, the rocky planet does not remain bound and hence does not 
reach the core. 

In the calculations presented so far, the mass and radius of the Jovian planet were taken 
to be those of Jupiter. However, many Hot Jupiters have infiated radii, with typical values 
Rp ~ l.2Rj for a mass Mp = 1 Mj. In addition, the radii of young g iant planets — befor e 



they contract and cool — are even larger than those of mature planets ( iBurrows et al.lll997l ). 
As a result, collisions of rocky planets with infiated gaseous planets (for the same mass) must 
be explored. To bracket the possibilities, the previous calculations were repeated for a Jovian 
planet with mass Mp = IMj, but with a larger radius Rp = 2Rj, leading to a lower mean 
density of only p ~ 0.1 g cm~^, ten times lower than the mean density of Jupiter. These 
results are shown in Figure [HI Compared to the previous case (with Rp = IRj), the impact 
energies are larger by an order of magnitude. Notice also that the curves do not extend to 
large impact parameters (bottom panel) because the rocky planets pass through the gaseous 
planet. 

Another way to illustrate the effects of inflated radii for the Jovian planets is to the plot 
the impact energy as a function of the radius Rp of the gas giant. These results are shown 
in Figure |9] for collisions with rocky planets that have starting speeds fo = 50 km s~^ and 
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Fig. 7. — Kinetic energy of the rocky planet as it reaches the core of the Jovian planet with 
radius Rp = Rj. Impact energy is shown as a function of initial velocity for fixed Xq = 1 
(top panel) and as a function of impact parameter at fixed impact speed Vq = 80 km s~^ 
(bottom panel). The four curves correspond to four masses for the rocky planet, for rrip = 
0.1, 1, 10, and 20 (from bottom to top). The Jovian planet has polytropic index n = 1, 
mass Mp = Mj, and dimensionless core radius = 1- Note that the curves do not extend 
beyond Xq ~ 0.9; for collisions with larger impact parameters, rocky planets pass through 
the gaseous planet and escape. 
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Fig. 8. — Kinetic energy of the rocky planet as it reaches the core of a young Jovian planet 
with large radius Rp = 2Rj. Impact energy is shown as a function of initial velocity for 
fixed Xo = 1 (top panel) and as a function of impact parameter at fixed impact speed f o = 
80 km s^^ (bottom panel). The four curves correspond to four masses for the rocky planet, 
for rrip = 0.1, 1, 10, and 20 (from bottom to top). The Jovian planet has polytropic 
index n = 1, mass Mp = Mj, and dimensionless core radius = 1- Note that the curves 
do not extend all the way to Xq = 1; for collisions with sufficiently large impact parameters, 
rocky planets pass through the gaseous planet and escape. 
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varying masses nip = 0.1, 1, 10, and 20 Mq. The kinetic energy Ef remaining when the 
rocky planet strikes the core of the Jovian planet is a smoothly increasing function of the 
radius of the larger body. To leading order, the four curves shown in Figure M have nearly 
the same shape, so that they are scaled by the mass of the rocky planet. The starting energy 
for the 1 planet is 7.5 xlO'^" erg; for collisions with Jovian planets with Rp = 2.1Rj, the 
impact energy is about two thirds of this starting value. Rocky planets with smaller masses 
lose a somewhat greater fraction of their initial energy, whereas sufficiently massive rocky 
planets {nip > lOM^) gain more energy by falling through the gravitational potential well of 
the giant planet than they lose through frictional dissipation. 

This difference in the energy dissipation between Jovian planets with different radii 
is remarkable. The planets with inflated radii tend to dissipate far less energy in their 
outer envelopes, and the percentage of the initial energy remaining when the incoming rocky 
planets reach the core can be as high as 70%. Note that this percentage represents the 
fraction of the rocky planet's total energy (kinetic and potential) at the surface of the Jovian 
planet. If we examine purely the kinetic energy, we find that the energy dissipation in the 
outer regions is low enough that the kinetic energy can actually increase during the inward 
trajectory, as it travels deeper into the potential well of the giant planet. 

The consequences of these collisions differ greatly for young Jovian planets with low 
densities and mature planets with higher densities. Mature planets are able to dissipate 
much more of the incoming energy near the surface, where heat is transported efficiently by 
radiation. In such cases, the energy of the impact is expected to radiate away quickly, and 
the structure of the giant planet should recover on short time scales. On the other hand, 
young planets with large radii tend to dissipate little energy in the radiative zones, and 
the majority of the impact energy is delivered to the core, where it becomes trapped and 
is available for long-term heating. The evolutionary stage of the Jovian planet is therefore 
an important factor in determining the long-term effects of these collisions. Young inflated 
planets will be far more susceptible to core enhancement and prolonged heating than mature 
planets. Similarly, Hot Jupiters will be more susceptible than giant planets farther from their 
stars. 



4. Tidal Disruption 

For purposes of flnding their orbits, the impinging rocky planets have thus far been 
treated as indestructible point particles. In this section we consider the possibility that 
these planets can be disrupted by tidal forces. If the incoming planet is destroyed, several 
outcomes are possible. Depending on the interior conditions of the giant planet, the rocky 
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Fig. 9. — Effects of the Jovian planet radius on collisions with rocky planets. Kinetic energy 
of the rocky planet as it reaches the core of the Jovian planet is shown as a function of 
Jovian planetary radius Rp. The four curves correspond to rocky planet masses rup = 0.1, 
1, 10, and 20 (from bottom to top). The Jovian planet has polytropic index n = 1, 
mass Mp = Mj, and dimensionless core radius = 1- For all cases shown, the initial speed 
(evaluated at the surface of the Jovian planet) is Vq = 50 km s~^ and the impact parameter 
xq = (head-on collisions). 
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planet remnants could either become accreted onto the core, or uniformly enrich the plane- 
tary interior with heavy elements. A full determination of the outcome is computationally 
expensive and is beyond the scope of this present work; in addition, it requires a greater un- 
derstanding of Jovian planet interiors (e.g., including convection) than is currently available. 
In this section we consider how tidal forces act on the incoming rocky planets and determine 
the likelihood of these bodies reaching the core (of the Jovian planet) relatively intact. These 
results allow us to estimate the maximum amount of mass and energy deposited at the core. 

At any point in the trajectory, the tidal force (per unit mass) exerted on the rocky 
planet can be approximated by the difference in the gravitational acceleration at the near 
and far surfaces, 

, ^ _ GMjr) GMjr + Ar) 
~ r2 (r + Ar)2 ' 

where Ar = 2rp (the rocky planet diameter). Note that the tidal forces arise due to the 
difference in the enclosed mass as well as the difference in the radial distance. Since the 
tidal forces are not spherically symmetric with respect to the rocky planet, this treatment 
introduces an approximation. Nonetheless, if the rocky planet is to withstand tidal forces, 
the tidal acceleration given by equation f l2^ must be less than the surface gravity of the 
rocky planet. Equating these quantities yields the critical condition for survival, i.e.. 



GM{r) GM{r + Ar) 



J.2 I J. _|_ 
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AF = ^ ' - Mr^Gp{r). (24) 



This expression can be simplified by expanding in {Ar)/r ^ r and keeping only the leading 
order terms. With these approximations, the expression ( l22l) for the tidal force takes the 
more familiar form 

4rpGM(r) 

Near the surface of the Jovian planet, the density is small enough so that the second term 
on the right hand side is small compared to the first. After dividing by rp, we obtain a rough 
criterion for the survival of the rocky planet. 

In other words, the incoming body will be tidally disrupted if its mean density is less than 
(about) four times the mean density of the giant planet. Note that equation (l25ll is valid 
near the surface of the Jovian planet, where the differential change in mass and the ratio 
(Ar)/r are relatively small. Here we consider rocky planets with radii rp ~ and giant 
planets with radii Rp ~ lO-R®, so that (Ar)/r ^ 0.2 at the surface and increases with 
decreasing radial position. As a result, the approximation scheme deteriorates as the rocky 
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planet moves inward. Here we use equation fl25l) to provide an order of magnitude estimate 
for the tidal stretching at the Hot Jupiter surface, but use the full form of equation fl22|) to 
assess the prospects for rocky planet survival (see below). 

Notice that the two sides of equation (!25|) are nearly equal for a pair of planets with the 
mean densities of Jupiter and Earth (where (p) ~ 1.3 g cm~^ and ~ 5.5 g cm~^, respectively). 
As a result, an Earth-like planet can be tidally disrupted by a giant planet like our Jupiter, 
but very well may survive if the radius of the Jovian planet is somewhat larger (so that its 
density is lower). Since both Hot Jupiters and young Jupiters (in any orbit) have inflated 
radii, incoming rocky planets are likely to survive tidal disruption, reach their centers, and 
increase the masses of their cores. 

The ratio of the tidal force to the surface gravity of the rocky planet can be calculated 
as a function of radial distance using equation fl25]) . Here we explore a variety of planetary 
properties. The Jovian planet mass was held constant at Mp = IMj and the radius varied 
over the range Rp = l — 2Rj, leading to average densities (p) ~ 0.1 — 1.0 g cm~'^. The results 
are shown in Figure [TD] for Rp = IRj and in Figure ITT] for Rp = 1.2Rj, where each curve 
represents a different rocky planet density. The latter value corresponds to a typical radius 



for o bserved Hot Jupiters, which are inflated relative to expectations (e.g., iLaughlin et al. 



20111 ). For even larger values of Rp (not shown), the densities of the Jovian planet are small 
enough that little tidal disruption t akes place. The exp ected range of (mean) densities for 
rocky planets is p ~ 2 — 10 g cm~^ (j Valencia et al.l 120071 ). with planet masses mp = 0.1-20 



Mq. Note that planets less massive than ~ IM^ are not expected to survive. Most of the 
calculations were performed for coreless polytropes of index n = 1; adding a core does not 
significantly alter the strength of the tidal forces (and hence these curves). 

Notice that the tidal force AF vanishes at radii r ~ Rp/2. The tidal force becomes 
negative for smaller radii, indicating that an orbiting planet will be first be stretched and 
subsequently compressed during its inward trajectory. In the center of the Jovian planet 
(where r — )■ 0), the compression reaches a maximum and the force given by equation 
formally approaches a well-defined limit 

limAF = -^Gp.rp. (26) 

Note that this is a (mathematically) formal result, because we are taking the limit r — )■ 
but the rocky planet has a finite size Vp. 

The area between the dashed lines (set where (AF)//p = ±1) in Figures fTOl and ITT] 
indicates the regions where the surface gravity of the rocky planet can withstand the tidal 
acceleration. If the ratio (AF)//p > 1, the outer layers of the rocky planet are stripped 
away; if (AF)//p < —1, the planet experiences devastating levels of compression. Here we 
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Fig. 10. — Ratio of tidal forces to rocky planet surface gravity /p as a function of radial 
distance (where /p = Grrip/rp^). The Jovian planet is taken to have no core, mass Mp = Mj, 
radius Rp = Rj, and n = 1. Results are shown for rocky planets with mass nip = 1 (top 
panel) and nip = 10 (bottom panel). In each panel, the curves correspond to different 
mean densities of the rocky planet, ranging from pp = 2 g cm~^ (top red curve) to pp = 10 g 
cm~^ (bottom black curve). The region between the dashed lines indicates radial locations 
where the rocky planet can survive. 
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Fig. 11. — Ratio of tidal forces to rocky planet surface gravity as a function of radial 
distance for an inflated gaseous planet. Here the Jovian planet is taken to have no core, 
mass Mp = Mj, radius Rp = 1.2Rj, and n = 1. Results are shown for rocky planets with 
mass nip = 1 (top panel) and m-p = 10 (bottom panel). In each panel, the curves 
correspond to different mean densities of the rocky planet, ranging from pp = 2 g cm~^ (top 
red curve) to Pp = 10 g cm~^ (bottom black curve). The region between the dashed lines 
indicates radial locations where the rocky planet can survive. 
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consider a rocky planet to survive whenever the condition |(AF)//p| < 1 is satisfied all the 
way to the core. All but the smallest planets are able to penetrate deep into the interior 
(at least 50 — 70% of the distance to the center) before experiencing severe tidal disruption; 
in many cases, the rocky planet reaches the core intact. As expected, the small planets 
{nip ~ O.IM®) must be very dense to survive the entire trajectory, and we exclude them 
in the remaining discussion. Here we focus on rocky planets with masses rup = 1, 10, and 
20 M^. For Jovian planets with the mean density of Jupiter (1.2 g cm~^), a rocky planet 
with rUp = IM^ must have density greater than pp ^ 5 — 6 g cm^'^ for survival; larger rocky 
planets with rrip = IOM0 can survive with somewhat lower density pp ~ 3 — 4 g cm~'^. 
On the other hand, for Jovian planets with sufficiently high densities (roughly, where the 
inequahty of equation [25] is violated), the tidal forces at the surface can exceed the rocky 
planet surface gravity, and such planets can be tidally disrupted before even entering the 
atmosphere. 

After a rocky planet is tidally disrupted, its fate remain unclear. The resulting rocky 
debris is unlikely to continue the journey toward the core, because smaller objects are more 
susceptible to frictional dissipation (and further tidal disruption). The rocky debris is likely 
to remain in the outer envelope of the giant planet and enrich the atmosphere, but its 
longer-term fate should be the subject of further study. The outcome will depend on, among 
other things, the opacity of the Jovian planet. For example, high opacities lead to convection, 
which could cause the rocky debris to be uniformly mixed into the gaseous planetary interior. 

In summary, we find that large rocky planets {rrip ~ lOM^) have a good chance of 
reaching the core intact; even planets with rUp ~ IM^ are able to survive in many cases, as 
long as the Jovian planet density is sufficiently low. As a result, these collisions provide a 
viable mechanism for allowing significant amounts of rocky material to be be accreted onto 
the cores of gaseous planets, especially those with infiated radii like young Jovian planets 
and much of the observed Hot Jupiter sample. 



5. Conclusion 



This paper explores collisions between Jovian planets and smaller terrestrial bodies in 
order to identify possible changes in the structure of the larger planets. This work focuses on 
the scenario where the Jovian planet has migrated inward and entered into a tight orbit, and 
subsequently experiences collisions with rocky planets that migrate later. However, many of 
the results apply to other parts of parameter space, e.g. , the possible accretion of cores by 



gaseous planets formed via gr avitational instability (e.g.. lBosslll997l : lBoleyll2009l : iBoley et al. 



20 111 : iHelled &: Shubertll2009l ). These collisions can affect the structure of the target giant 
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planets in three ways: [a] The metalhcity is increased, [b] The core size can be increased, 
and [c] A substantial amount of energy is imparted to the giant planet. For Hot Jupiters, 
these impacts could thus help explain the observed anomalous radius distribution. 

Rocky planets are efficiently captured by giant planets and thereby increase the metal- 
hcity of the larger body. For sufficiently high impact speeds vq and large impact parameters 
Xq > x^, the projectile planet can pass through the giant planet and emerge out the other 
side. For nearly head-on collisions, the initial velocity required to to escape capture is 
much higher (v-npr^ ~ 1 0^ km s~^) than the expected impact speeds f o ~ 40 — 150 km s~^ 



( iKetchum et al.ll2011bl ). With these latter speeds, incoming planets can escape only for large 
impact parameters xo^O.SRp, where the densities are low. In other words, penetration is 
possible when the incoming rocky planet skims the surface of the Jovian target. As a result, 
most rocky planets lose enough kinetic energy through their initial passage that they remain 
gravitationally bound and are captured (see Figure |3l an analytic estimate of the required 
penetration speeds is given in Appendix |A]). Nonetheless, the capture cross section of the 
giant planet can be reduced; for typical collision speeds vq ~ 100 km s~^, the cross section 
is smaller than the geometrical area by ~ 30%. 

These collisions provide a plausible mechanism for explaining the massive cores that 
have been inferred for observed Hot Jupiters. Although metallicity increases are essentially 
automatic through such collisions, core masses can only increase if the rocky planets avoid 
tidal disruption and survive to reach the central regions. The expected parameter space 
includes cases where the impinging planets survive and cases where they are completely 
disrupted. As a benchmark, rocky planets with the mass and radius of Earth live near the 
threshold for tidal destruction when they collide with giant planets having the mass and 
radius of Jupiter. However, sufficiently dense projectile planets or inflated giant (target) 
planets allow for the rocky planets to survive all the way to the central regions (see Figures 
Hn] and [TT]) . This flnding is relevant because such collisions are expected to take place early 
in the system's evolution when the giant planets are larger, and also because Hot Jupiters 
are subject to additional heating sources and are thus larger in radius for a given mass. 
As noted above, this process also provides a mechanism through which planets formed via 
gravitational instability can acquire cores. 

For collisions to signiflcantly alter the metallicity or the core mass of a giant planet, it 
must accrete rocky bodies that represent total masses AM ~ 10 — lOOM®. The accumulation 
of such mass increments could take place through collisions with many smaller rocky bodies 
or via fewer larger bodies. Large rocky bodies, with mass rrip > lOM^, are more effective 
for several reasons: First, larger rocky bodies are subject to less eccentricity damping during 
migration, and are thus more likely to collide with the giant planets. More speciflcally. 
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eccentricity damping will be small enough to allow collisions when th e rocky planets are 



large enough to partially clear gaps (lArtymowicz 



993 



Kiev et al. 



2004h. which occurs when 



their Hill sphere exceeds the disk scale height (iCrida et al.l 120081 : iPapaloizou &: Terquem 
20061 ). These considerations indicate that relatively large rocky planets (with mass rrip 
~ 10 — 20M^) are favored in order to achieve partial gap-clearing, reduced eccentricity 
damping, and ultimately collisions. In addition, collisions with larger rocky bodies allow for 
the required mass increments to be realized through a lower number of encounters. Finally, 
larger rocky bodies are more likely to survive tidal disruption and reach the central regions 
(see Figures [10] and [TTl) . On the other hand, if the circumstellar disk h as high levels of 



turbu lence, small rocky planets are subject to stochastic migration (e.g., lAdams &: Bloch 
20091 ). and c an still experience enc ounters with gas giants. Such encounters often result 
in collisions (IKetchum et al.ll2011b( ) when they take place at small semimajor axes (a<0.1 
AU); for larger a > 1 AU, most encounters resul t in ejection of the sma ller rocky planet, with 
collisions taking place only ~ 10% of the time (IKetchum et al.ll2011al ). 



Collisions provide a substantial amount of energy to the target giant planet. For rocky 
planets with mass rup = lOM^ and initial speed Vq = 60 km s~^, the kinetic energy trans- 
ferred to the giant planet is AE ~ 10^^ erg. If, for example, this energy is radiated over a time 
scale of 1 Gyr, the associated increment in planetary luminosity would be ALp 4 xlO ^^ 



WjWhich is large enough to affect the internal structure of the planet (iBodenheimer et al. 



20031 ) . However, in order for these collisions to provide an effective (long-term) heating 
source, most of the energy must be deposited deep in the Hot Jupiter interior, in regions of 
high opacity. For a given rocky planet mass, the amount of energy deposited depends on 
the density of the Jovian planet. For a Jovian planet with the mass and radius of Jupiter, 
most of the kinetic energy is dissipated in the outer layers of the planet. At the other end 
of the possible range, for a Jovian planet with the mass of Jupiter and an inflated radius 
Rp = 2Rj, little energy is dissipated in the outer layers and a sizable fraction remains upon 
impact with the core (see Figures O and [HI and Appendix [Bl). 

An important unresolved issue is the question of how quickly the collision-induced energy 
is transferred out of the planet. Collisions with large terrestrial planets deposit enough energy 
in the central regions of the Jovian planets to potentially affect their structure. This energy 
could even play a role in explaining the radius anomalies of the transiting exoplanets if 
the time scales for energy loss are sufficiently long. In future work, interior models of Hot 
Jupiters should be modified to include this possible heat source. 

As shown herein, incoming rocky b odies are so r netim es tid ally destroyed a fter entering 
the surface of the giant planet (see also lAnic et al.l (120071 ) and iLi et al.l (120101 ) for detailed 



simulations of specific cases). This work indicates that most planets heavier than nip ~ lOMg 
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are able to withstand tidal disruption (see Figures [TO] and [TT]) . However, this work focuses 
on frictional forces and tidal disruption, but does not consider the possible destruction of 
incoming rocky planets through heating and ablation. Through both ablation and tides in- 
coming planets continuously lose mass throughout their descent, and the relative importance 
these effects should be studied in future work. If ablation is strong enough, the rocky planets 
could completely disintegrate before reaching the core, even if they can survive tidal disrup- 
tion. Ablation will thus increase the overall metallicity of the Jovian planet atmosphere and 
reduce the amount of kinetic energy and mass delivered to the core. Since the calculations 
of this paper do not include ablation, the quantities of heat and heavy elements deposited 
deep in the Jovian planet interior could be lower than estimated herein. 

Although ablation effec ts are important for smaller rocky bodies such as planetesimals 



Anic et al. 


2007; 


Li et al. 


2010) 



planets with masses rrip = 1 - 10 i\ 
incoming projectile planet is expected to be destroyed only after it encounters a gas mass 



larger than its own mass (IKorycansky fc Zahnld 120051 ). The mass in gas Am encountered 



on an inward trajectory is roughly given by Am = UpppRp, so that survival requires a 
minimum size > SppRp/^App) ~ 1 — 2R^. As a result, the larger rocky planets of interest 
here are likely to survive ablation, but smaller bodies are susceptible to destruction. 

The results of this paper show that the mean density of the giant planet plays a defining 
role in determining the consequences of these collisions. The average density helps determine 
whether or not incoming rocky planets can survive tidal disruption and affects the radial 
layers where kinetic energy is deposited. Since giant planets are born with enlarged radii, 
this density dependence translates into a dependence on the age of the planet. A rocky 
planet that collides with a young giant planet (with low density) can have dramatically 
different consequences than a rocky planet that collides with a mature giant planet (with 
higher density). The environment in which the giant planet resides plays an additional role. 
Since Hot Jupiters often remain significantly inflated throughout their evolution, this subset 
of extrasolar planets should be even more sensitive to these collisions. This lends promising 
support for the proposal that planetary impacts could help explain, in part, the anomalous 
radius distribution of the observed Hot Jupiters. 
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A. Analytic Estimate for the Penetration Velocity 

In this Appendix, we derive an analytic estimate for the penetration velocity Vpen, the 
speed required for an impinging rocky planet to pass through a Hot Jupiter and thereby 
avoid capture. As shown below, the required velocity is much larger than the gravitational 
escape speed fesc from the planet. As a result, gravity can be neglected to leading order, 
and the equation of motion simplifies to the form 

^ = -aPv'v = -arv{xx + ijy) , (Al) 

where all of the variables are dimensionless (defined in the text). By definition, the initial 
velocity lies in the y direction and the initial spatial position is given by yo = —L, where 
L > 0. Further, since the required speeds are large, the path of the rocky planet through 
the Hot Jupiter can be approximated as a straight line (for cases of interest where the rocky 
planet is not captured). As a result, we can set x ~ and v = y. With these simplifications, 
equation (lAll) becomes 

nil 

- = -/3^;^ /3 = a(r), (A2) 

where we have made the approximation /" = (/"). This differential equation is easily 
integrated. Applying the boundary condition v = Vq at time r = 0, we obtain the result 

'=d^ = TT^o- ^^^^ 

Integrating equation (1A3P for ?/(r) yields the solution 

y-yo = r'H^ + voM- (A4) 
We can combine equation (1 A3 P and flA4p to find the speed f as a function of position y, 

v = voe-^^y-y°l (A5) 

In order for the rocky planet to penetrate the Jovian planet, the speed must exceed the 
escape speed (f > i^esc) at the location y = +L = —yQ. The minimum initial speed fpen 
required for penetration is thus given by the condition 

^esc ^pen 6 , Or t'pen ^esc 6 (A6) 

Note that the escape speed is given by t>esc = (2/Uo/^o)^^^ in terms of dimensionless variables. 

Next we find the dependence of the penetration speed fpen on the impact parameter 
xq. Here we estimate the average density by approximating the density profile of the Jovian 
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planet as a Gaussian function, where the peak is a function of the impact parameter. For 
head on colhsions, with xq = 0, we assume (/") ~ 1/2. The constant /3 is thus given by 

P = ^a exp [-7 Xq] , (A7) 

where we have introduced an unspecified parameter 7, which is expected to be of order unity. 
We find that reasonable agreement with numerical results can be obtained using 7 1/4 
(see Figure [3]) . Equations (1A6P and (]A7|) together yield a simple relation between the initial 
velocity and the impact parameter. In the limiting case of grazing collisions, where xq = ^0 
and 1/0 = 0, the penetration speed is equal to the escape velocity fpen = Vesc- On the other 
hand, for head-on collisions where Xq = and z/q = —^q, the penetration speed attains its 
maximum value fpen ~ t^esc exp[a^o]- Note that even though the parameters a and ^0 are 
"of order unity" , their typical values are ~ 3, so that the exponential factor is large and 

^pen ^esc- 

The above derivation assumes that (f^) ^ 1/2. For the particular case of an n = 1 
polytrope with no central core (^c — ?■ 0), we can directly test this approximation: 

1 n sin 1 
{n = {f) = - — ^rfe = -Si(vr) , (A8) 

where Si(x) is the sine integral (lAbramowitz fc StegunI 119701 ). which has the known value 
Si(7r) 1.85194, so that (/") k, 0.5895. Thus, the approximation is reasonable. 



B. Analytic Estimate for Energy Dissipation 

In this Appendix, we provide an analytic estimate for the energy Ef remaining when an 
incoming rocky planet reaches the core of the Jovian planet. The incoming planet has initial 
kinetic energy £^0 = '^^o/'^^ where the speed vq is evaluated when the rocky planet reaches the 
surface of the Jovian planet. For the particular case of direct collisions, with zero impact 
parameter, the incoming planets travel along a radial path from the gas giant surface to the 
core. The amount of energy W dissipated along this path is given by 

W = a r'rv'd^, (Bl) 

where we have chosen signs so that W > 0, and where all quantities are dimensionless (see 
text). Using equation flA5|) to estimate the speed as a function of position, the work W can 
be expressed as 

W = avo" [ "/"e'''^^-^") , (B2) 



- 33 - 



where (5 = a{f^). If we make the simphfying approximation /" ^ 1/2, the expression 
becomes ^ 

W = -avo^ [ "e^^-^") d^, (B3) 

which can be integrated to yield 

W = ^wo' [1 - e"(«=-«°)] = ^0 [1 - e"(^^"«°)] . (B4) 

Conservation of energy imphes Ej = Eq — W, so the remaining energy when these planets 
reach the core is given by 

Ef = e'^^^^-^"^ Eo . (B5) 

Consider a giant planet with polytropic index n = 1 and core radius = 1- The final 
energies Ef predicted by equation ( ]B5l) for incoming rocky planets with masses rup = 1, 
10, and 20 are given by Ef/Eo ~ 0.03, 0.20, and 0.30, in rough agreement with the 
numerically determined results. 
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